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Abstract. We establish general product formulas for the solutions of non- 
autonomous abstract Cauchy problems. The main technical tools are evo- 
lution semigroups allowing the direct application of existing results on au- 
tonomous problems. The results obtained are illustrated by the example of 
an autonomous diffusion equation perturbed with time dependent potential. 
We also prove convergence rates for the sequential splitting applied to this 
problem. 



Operator splitting procedures are used to solve ordinary and partial differential 
equations numerically. They can be considered as certain finite difference meth- 
ods which simplify or even make the numerical treatment of differential equations 
possible. The idea behind these procedures is the following. In many situations, 
a certain physical phenomenon can be considered as the combined effect of several 
processes. Hence the behavior of a physical quantity is described by a partial dif- 
ferential equation in which the time derivative depends on the sum of operators 
corresponding to the different processes. These operators usually are of different 
nature and for each sub-problem corresponding to each operator there might be an 
effective numerical method providing fast and accurate solutions. For the sum of 
these operators, however, it is not always possible to find an adequate and effective 
method. Hence, the idea of operator splitting procedures means that instead of the 
sum we treat the operators separately and the solution of the original problem is 
then to be recovered from the numerical solutions of these sub-problems. We refer 
to the recent monographs by Farago and Havasi [TT] or Holdcn et al. [15] for a 
detailed introduction to the theory and applications of operator splitting methods. 

There was enormous progress in recent years in the theoretical investigation of 
operator splitting procedures. Especially, ordinary differential equations and au- 
tonomous linear evolution equations have been treated thoroughly, sec also Batkai, 
Csomos and Nickel [5] and the subsection below for a (certainly not complete) list 
of references. 

The aim of the present paper is to investigate the above described splitting 
method for non-autonomous evolution equations of the form 
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on some Banach space X. Our particular goal is to emphasize that non- autonomous 
evolution equations can often be rewritten as an autonomous abstract Cauchy prob- 
lem by means of an appropriate choice for the state-space. Thus, by making use 
of so-called evolution semigroups, it is possible to apply existing results for au- 
tonomous problems. 

First we summarize the necessary background on operator splitting for abstract 
Cauchy problems, i.e., operator splitting in the framework of strongly continuous 
operator semigroups. The key ingredient here is Chernoff 's Theorcm ll.il Then non- 
autonomous evolution equations and evolution semigroups arc surveyed, providing 
the main technical tools for the succeeding sections. 

A product representation is presented in Section [3J while operator splitting — 
strictly in the sense above — is considered in Section [3] To keep our presentation 
short, we mainly restrict ourselves to the case of the so-called sequential splitting, 
but in Section|4]we show how higher order splitting methods can be treated with es- 
sentially no difference. In that section, we also prove the convergence of the splitting 
methods when combined with spatial "discretization," and make a quick outlook 
on the positivity of evolution families. Finally, as an illustration of the developed 
tools, we apply them to a diffusion equation with time dependent potential. More- 
over, by semigroup methods, using results of Jahnke and Lubich |18j . and Hansen 
and Ostcrmann |13l 114] . we obtain estimates on the order of the convergence. 

A word on notation: For a family of operators Uq, U\, . . . , U n ~\ 6 5£{X), we 
denote the ( "time-ordered" ) product of these operators by 
n-l o 

l[U p :=U n - 1 U n -2---U 1 Uo and J[ U p := U Q U X ■ ■ ■ U n ^U n ^. 

p—Q 

Operator splitting for autonomous problems. In this section, we recollect the 
main notions and results of operator splitting for autonomous equations. Consider 
the following abstract Cauchy problem on a given Banach space X: 



(ACP) 




= (A + B)u(t), t>0, 
= x€X, 



where the operators A, B, and the closure C := A + B are supposed to be gen- 
erators of strongly continuous semigroups T, S, and U, respectively. Our general 
reference on strongly continuous operator semigroups is the monograph Engel and 
Nagel 0. 

As mentioned in the introduction, operator splitting means that we try to recover 
the solution semigroup U using the semigroups T and S. As for splitting procedures 
we mention the most frequently used ones (for more details, see Batkai, Csomos 
and Nickel Section 2.2]): 

• The sequential splitting, classically the Lie- Trotter product formula, is given 

by 

u?(t) := [S(t/n)T(t/n)] n Xl 

• the Strang splitting is given by 

uSt(t) : = [T(t/2n)S(t/n)T(t/2n)] n x, 

• and — for a fixed parameter O £ (0, 1) — the weighted splitting is 

<(t) := [QS(t/n)T{t/n) + (1 - Q)T{t/n)S(t/n)] n x 

with n € N. In case O = h, it is also called symmetrically weighted splitting. The 
convergence of these procedures is usually ensured by the following classical result. 
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Theorem 1.1 (Chernoff [BJ, or see Engel and Nagel [5J Sec. III.]). Let C be a linear 
operator in the Banach space X and assume that F : R + — > _Sf (X) is a (strongly) 
continuous function with F(0) = I and 

\\(F(t)) k \\ < Mc kult for allt>0 andkeN (stability). 

Suppose that there is a dense subspace D, with (A — C)D being also dense for some 
(large) A > 0. If for every x € D the limit 

lim — = Cx (consistency) 

exists, then C is the generator of a Co-semigroup U , the set D is a core for the 
generator C , and we have 

lim (F( — ))™x = U(t)x (convergence). 

Note that if the closure of C is already known to be a generator, as it is the case in 
problems motivated by numerical analysis, then the range condition is automatically 
satisfied. 

The operator family F is sometimes called a finite difference method. Clearly, 
the above mentioned splitting procedures have this form. For example, for the 
sequential splitting we take 

F sq (h) = S{h)T(h). 

It is important to note that Chernoff 's Theorem does not yield anything a priori 
about the rate of convergence. The finite difference method F is said to be of order 
p > 0, if for x from a suitably large subset of X there is C > such that for all 
t £ [0, to] we have 

\\FG) n x-U(t)x\\<&, 
or, as in many special cases, equivalently, 

\\F{h)x-U{h)x\\ < C'h p+1 . 

The equivalence holds in special cases where it is possible to ensure the invariance 
of the above mentioned large subset D of X (for more details we refer to the Lax 
equivalence theorem which states that the above two definitions are equivalent for 
a finite different method if and only if the method is stable). 

Different splitting procedures were introduced to increase the order of conver- 
gence. In the finite dimensional setting, it is well known that the sequential splitting 
is of first order, the Strang and the weighted splitting with = i are of second 
order. Moreover, the weighted splitting allows also the use of parallel computing. 

In the infinite dimensional case, however, no similar general statement can be 
made without additional assumptions. There has been intense research in this 
direction, and we mention the works by Bj0rhus [4], Cachia and Zagrebnov [5], 
Farago and Havasi [THI, Hansen and Ostermann [13], Ichinose et al. [17], Jahnke 
and Lubich [T5] or Neidhardt and Zagrebnov [27] . 

To obtain error estimates later for diffusion problems, we apply a result by 
Jahnkc-Lubich, Hanscn-Ostcrmann, which relies on commutator bounds. For sim- 
plicity, we mention here only the special case used later. 

Theorem 1.2 (Jahnke and Lubich [181 Theorem 2.1], Hansen and Ostermann 
[131 Theorem 2.3]). Suppose that A generates a strongly continuous contraction 
semigroup e tA in the Banach space X and that B € ^f(X) such that there exists 
an a > such that 

(1) \\[A,B]v\\ = \\(AB - BA)v\\ < c||(-A)%|| 
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for all v £ D C D((—A) a ) (where D is some dense subspace of D((—A) a ) invariant 
under c t( - A+B 1). Then one has first order convergence for the sequential and Strang 
splittings, i.e., 



e t(A+B) v 



< 



< 



n 

ce 



(-A) a v 
{-A) a v 



Non-autonomous evolution equations and evolution semigroups. In this 
section we summarize the main results and definitions on non- autonomous evo- 
lution equations and evolution semigroups needed for our later exposition. For a 
detailed account and bibliographic references see, e.g., the survey by Schnaubelt in 
[SI Section VI. 9.]. Consider now the non-autonomous evolution equation 



(NCP S , X ) 



' ±u(t) = A(t)u(t), t>se 

u(s) = x G X, 



where X is a Banach space, (A(t),D(A(t))) is a family of (usually unbounded) 
linear operators on X. 

Definition 1.3. A continuous function u : [s, oo) — > X is called a (classical) 

solution of (NCP^) if u e C 1 ([s, oo); A"), u{t) g D{A{t)) for all t > s, u{s) = x, 
and jtu{t) = A(t)u(t) for t > s. 

We use the following slight modification of Kellcrmann's definition Definition 
1.1] for the well-posedness of the non-autonomous Cauchy problem (NCP). 

Definition 1.4 (Well-posedness). For a family (A(t), D(A(t))) R of linear op- 
erators on the Banach space X the non-autonomous Cauchy problem (NCP) is 
called well-posed (with regularity subspaces (Y s ) s6 i{ and exponentially bounded 
solutions) if the following are true, 
(i) (Existence) For all s € K the subspace 



Y. 



there exists a classical solution for 



r (NCP)^} c D(A(s)) 



is dense in X. 

(ii) (Uniqueness) For every y €Y S the solution u s (-, y) is unique. 

(iii) (Continuous dependence) The solution depends continuously on s and y, 
i.e., if s n — > s € K, y n —ty£Y s with y n G Y Sn , then we have 

\\u Sn (t,y n ) - u s (t,y)\\ -> 
uniformly for t in compact subsets of R, where 



u r (t,y) 



u r (t,y) if r < t, 
y if r > t. 



(iv) (Exponential boundedness) There exist constants M > 1 and w£l such 
that 

\\u s (t,y)\\ < Afe^*— 5||y|| 
for all i/Efs and t > s. 

As in the autonomous case, the operator family solving a non-autonomous Cauchy 
problem enjoys certain algebraic properties. 

Definition 1.5 (Evolution family). A family U = (U (t, s))t> s of linear, bounded 
operators on a Banach space X is called an (exponentially bounded) evolution 
family if 
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(i) U(t, r)U{r, s) = U{t, s), U{t, t) = I holds for all i > r > s G R, 

(ii) the mapping (t, s) H> U (t, s) is strongly continuous, 

(iii) \\U(t, s)\\ < Mc^ 1 "^ for some M > l,u G R and all t > s G M. 

In general, however, and in contrast to the behavior of Co-semigroups (i.e., the 
autonomous case), the algebraic properties of an evolution family do not imply 
any differentiability on a dense subspace. So we have to add some differentiability 
assumptions in order to solve a non-autonomous Cauchy problem by an evolution 
family. 

Definition 1.6. An evolution family U = (U(t, s)) t > s is called evolution family 

solving (NCP) if for every s G R the regularity subspace 

Y a ■= [y G X : [s, oo) 3 t h-> U (t, s)y solves (NCP) s y | 
is dense in X. 

The well-poscdness of (NCP) can now be characterized by the existence of a 
solving evolution family. 

Proposition 1.7 (Nickel [^HJ Proposition 2.5]). Let X be a Banach space, and 
assume that (A(t), D(A(t))) R is a family of linear operators on X and consider the 
non- autonomous Cauchy problem (NCP). The following assertions are equivalent, 
(i) The non- autonomous Cauchy problem (NCP) is well-posed, 
(ii) There exists a unique evolution family (U(t, s))t> s solving (NCP). 

To every evolution family we can associate Co-semigroups on JT-valued function 
spaces. These semigroups, which determine the behavior of the evolution family 
completely, are called evolution semigroups. Consider the Banach space 

BUC (R; X) = {/ : R — > X : / is bounded and uniformly continuous}, 

normcd by 

ll/H :=sup||/(t)||, /GBUC(R;X); 

tern 

or any closed subspace of it that is invariant under the right translation semigroup 
7Z defined by 

(K(t)f)(s) := f(s - t) for / G BUC (R; X) and s G M, t > 0. 

In the following X will denote such a closed subspace; we shall typically take X = 
Co(R; X), the space of continuous functions vanishing at infinity. 

It is easy to check that the following definition yields a strongly continuous 
semigroup. 

Definition 1.8. For an evolution family U — (U (t, s)) t > s we define the correspond- 
ing evolution semigroup T on the space X by 

(T(t)f)(s) ~U(s,s-t)f(s-t) 

for / G X, s G R and t > 0. We denote its infinitesimal generator by (Q,D(Q)). 

With the above notation, the evolution semigroup operators can be written as 

T(t)f = U(;--t)K(t)f. 

We can recover the evolution family from the evolution semigroup by choosing a 
function / G X with /(s) = x. Then we obtain 

(2) U(t >S )x = (K(s-t)T(t-s)f)(s) 

for every s G R and t > s. 
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The generator of the right translation semigroup is essentially the differentiation 
— 4- with domain 

as 

£>(-&) := X, := {/ G C^R;*) :f,f'£ X). 

For a family (A(t),D(A(t))) t g of unbounded operators on X we consider the 
corresponding multiplication operator (A(-), Z?(A(-))) on the space X with domain 

L>(A(-)) := {/ € X : f(s) £ U(i(s))Vs£ R, and [s h> A(s)/(s)] £ A'}, 
and defined by 

(A(-)f)(s) := A(s)f(s) for all s £ R. 
Now we characterize well-posedness for non-autonomous Cauchy problems. 

Theorem 1.9 (Nickel Theorem 2.9]). Given a Banach space X , and a family of 
linear operators (A(t), D(A(t))) teR on X. The following assertions are equivalent. 

(i) The non- autonomous Cauchy problem (NCP) for the family (A(i)) tS R is well- 
posed (with exponentially bounded solutions), 
(ii) There exists a unique evolution semigroup T with generator (Q, D(Q)) and an 
invariant core T> C X\ PI D(Q) such that 

Qf + f' = A{-)f 

for all f £ V. 

Conditions implying well-posedness are generally divided into assumptions of 
"parabolic" and of "hyperbolic" type. Roughly speaking, the main difference be- 
tween these two types is that in the parabolic case we assume all A(t) being genera- 
tors of analytic semigroups, while in the hyperbolic case we assume the stability for 
certain products instead. In both cases one has to add some continuity assumption 
on the mapping t i— > A(t). We mention only a typical and quite simple version for 
each type. 

Assumption 1.10 (Parabolic case). 

(PI) The domain D := D(A(t)) is dense in X and is independent of t £ R. 

(P2) For each ( 6 R the operator A(t) is the generator of an analytic semigroup 

e -A(t)_ For all t g Kj the re solvcnt i?(A, A(t)) exists for all A £ C with 5RA > 

and there is a constant M > 1 such that 

\\R(X,A(t))\\<j^i 

for 5ft A > 0, t £ R. The semigroups e' A W satisfy \\e sA ^\\ < Me us for absolute 
constants u> < and M > 1. 
(P3) There exist constants L > and < a < 1 such that 

\\(A(t) -^(s))A(0)" 1 !| <L\t~s\ a for all £, s e M. 

Assumption 1.11 (Hyperbolic case). 

(HI) The family (A(t)) t ^M is stable, i.e., all operators A(t) are generators of Co- 
semigroups and there exist constants M > 1 and w £ R such that 

(uj, oo) C p{A(t)) for all * £ M 

and 

k 

|JJ R(X,A(tj)) <A-/(A-cj)~ fc forallA>u; 

3=1 

and every finite sequence — oo < t\ < t-2 < • • • < t^ < oo, fceN. 
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(H2) There exists a densely embedded subspace Y X, which is a core for every 
A(t) such that the family of the parts (A|y(i)) t6 R in Y is a stable family on 
the space Y. 

(H3) The mapping R3ii-> A(t) £ C(Y,X) is uniformly continuous. 

Remark 1.12. Since the classical papers of Evans [S], Howland [16], and Neid- 
hardt gU [13 [55] , evolution semigroups have been intensively used to study non- 
autonomous evolution equations. Here, various results on well-poscdness as well as 
qualitative behavior of these equations were obtained. For a quite comprehensive 
overview and a long list of different variants we refer, e.g., to Nagel and Nickel 
P). Neidhardt and Zagrebnov gS], Nickel [Ml [3D], Nickel and Schnaubelt [3T] and 
Schnaubelt [33] . The recent article Neidhardt and Zagrebnov [35] focuses on (quite 
general) assumptions of the "hyperbolic" type and obtains well-posedness results 
- in a general sense - for non-autonomous evolution equations by properly defin- 
ing and analyzing the "sum" — ^ + A(-) yielding the generator of the associated 
evolution semigroup. In contrast to that approach, in our paper wc simply as- 
sume well-posedness of our evolution equation under any appropriate (parabolic or 
hyperbolic) condition. Therefore, the solving evolution family, the corresponding 
evolution semigroup, and its generator are well defined by assumption. Our main 
interest is then, how these solutions can be approximated (numerically) by splitting 
procedures. 

2. A PRODUCT FORMULA 

In this section we present a product formula for the solutions of the non-autono- 
mous Cauchy problem (|NCP[) . In the case B(t) = 0, this formula essentially goes 
back to Kato [T5]. This splitting-type formula is especially useful if for every time 
r£M we are able to solve effectively the autonomous Cauchy problems 

(Eq. 1) £«(t) = A(r)u(t) 

(Eq. 2) ^v(t) = B(r)v(t) 

with appropriate initial conditions. This is usually the case if the operators A(-) 
and -B(-) are partial differential operators with time dependent coefficients or time 
dependent multiplication operators. Formally, this means that we assume that the 
operators A(r) and B{r) generate strongly continuous operator semigroups, which 
we denote by using the exponential notation as e' A ^ and e' B ^ r \ respectively. We 
devote this section to the simplest product formula arising from the sequential 
splitting. 

Suppose we want to determine the solution of (NCP) at time t + s > and hence 
take the time-step r = t/n. We start with the known initial value u sq (s) = x, then 
solve the first (Eq. 1) equation on the time interval [s, s + r] taking r = s. Then 
we take the result ( s + T ) a s the initial value for the second equation (Eq. 2) 
which we solve again on [s, s + t]. With this new result u sq (s + t) := itj (s + t) as 
initial value for (Eq. 1) we restart the procedure and iterate it n times. Formally: 

A(s + {k~l)T)u[ k) {t), te (a + (ft - l)r, s + At] , 
w sq (s + (fc - l)r), 

B(s + (fc-l)T)4 fc) W: *G (s + (fc - I)r, s + fcr] , 
u[ k) {s + kT), 

4 k \ s + k T ), 



u[ k \s + (k - l)r) = 

u [k \s + (k l)r) = 
u sq (s + kT) := 
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with k = 1, 2, . . . , n. Using that for r £ [0, r], 

u[ k] (s + (k - 1)t + r) = e rA{s+{k - 1)T) u sq (s + (k - 1)t), 

and that 

4 fe) (s + (k - l)r + r) =e rB ^ k - 1 ^u[ k) {s + kr) 

=c rB( S +(fe-l)r) e rA( S +(fc-l)r) uSq ^ + ^ _ ^ 

wc see by a simple induction argument that the split solution u sq (s + /ct), obtained 
by applying the sequential splitting procedure, can be written as 

fc-i 

(3) u sq (s + fcr) = Y[ c TB{ - s+pT h TA{s+pT) x for fc e N, kr < t, and x £ X. 

In what follows, we study the convergence of this expression. 
Assumption 2.1. Suppose that 

a) the non-autonomous Cauchy problem corresponding to the operators (A(-) + 
£?(•)) is well-posed, 

b) (Stability) the operators Air) and Blr) are generators of Co-semigroups e' A ( r \ 
Q B ( r ) of type (M, uj) (M > 1 and w £ 1) on the Banach space X and, therefore, 

(w, oo) C p(A(r)) n p(B(r)) for all r£l. 

Moreover, let 



sup 1 1 | | ^e™ ™ 'e 



< Me"', and 



c) (Continuity) the maps 

t ^ R(X, A(t))x, fn> R(X,B(t))x 
are continuous for all A > cj and a; G X. 

We denote the evolution family solving (|NCP|) by W and the corresponding 
evolution semigroup, generated by the closure C of C := — 4- + A(-) + B(-), by W. 

As we shall see in a moment, Assumption 12.11 yields that the multiplication 
operators A(-), B(-) with appropriate domain generate strongly continuous multi- 
plication semigroups on Co(R; X) (for more on this matter we refer to Engel and 
Nagel [El Sec. III.4.13] and Graser QJ] ). 

Theorem 2.2. Under Assumvtion \2. 1\ one has the convergence 

n-l 

(4) W(t,s)x= lim TT (e^'+^e^'+^Ax 

p=0 



for all x £ X, locally uniformly in s,t with s <t. 

Proof. The main idea of the proof is analogous to the one in Nickel [30l Proposition 
3.2]. Consider the semigroups e' A ^ and c' B ^ for given r £ R. By the uniform 
growth assumption in 12. li b) on the semigroups, for fixed t > the function r t— > 
e tA ^f(r) vanishes at infinity whenever / has this property. We also have that the 
function r i-> e tA ^ is strongly continuous. Indeed, by the Trotter-Kato Theorem 
(see Engel and Nagel [8j Thm. III.4.8]) we even obtain that K+ x R 3 (t, r) i-> e M(r) 
is strongly continuous. All these reasonings are, of course, true if Air) is replaced 
by B(r). Let now / £ BUC(R;X). Then r h-> e' A M/(0 is continuous, too. We 
have therefore shown that the multiplication semigroups e tA ^ and e tB ('\ generated 
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by the multiplication operators A(-) and £>(•), both act on the space X = Co(K; X), 
see also Graser [T3]. It can be seen by induction that 

(*(A)e**<V*->) w /(0 = U(e^-^e^-^)n(t)f(.). 

p—n 

The stability assumption I2.1l b) immediately implies the stability for the finite 
difference method F(h) := TZ(h)c hB ^c hA ^'\ Consistency is standard to check: take 
/ G A?i n D(A(-)) n D(B{-)). Then we can write 



lim — - — — — = lim 

hlO h hlQ 



, C 



hA( 



K{h)e hB (->- i- + 



= A(-)f + B(-)f-f. 

By our well-posedness assumptions, the closure of the operator C = — -r- + !?(•) + 
A(-) generates a strongly continuous semigroup on X, hence the set (A — C)D(C) 
is dense in X. By the stability assumption we can apply Chernoff's Theorem 11.11 
with the three operators — 4-, A(-), B(-), and obtain that the evolution semigroup 
generated by C is given by 

W(t)/ = lim n( e ^C--£)e^--£>)/(.-i). 

p—n 

The above limit is to be understood in the topology of X, that is, in the uniform 
topology. By using this, and by applying the formula (|2|) from the previous section, 
we can recover the evolution family from the evolution semigroup and arrive at the 
formula 

W(t t s)x= lim rf(e¥^('-^) e ¥^-^)) Ii 

p—n 

from which the assertion follows. □ 

Remark 2.3. In the proof of Theorem [22] we have used that the semigroups e' A ^ 
and e' B ^ map Co(K;X) into itself. If c' A ^ and c' B ^ are uniformly strongly 
continuous in r G K, then one could also work on the space X = BUC (R; X). 

Remark 2.4. The stability condition b) is automatically satisfied, if A(t) and B(t) 
are generators of quasi-contractive semigroups with uniform exponential bound w 
for all t. 

Remark 2.5. In Vuillcrmot et al. [351121]: the authors prove the representation 
formula (0]) where A{t) and B(t) are generators of contraction semigroups, the 
family A(-) satisfies a version of the so-called parabolic condition and the family 
£?(•) is a small perturbation. Theorem 12.21 can be seen as a generalization of this 
result and can be applied not only in a larger class of parabolic problems but also 
in the hyperbolic case. In |35| Vuillermot proves a Chernoff-type approximation 
theorem for time-dependent operator families. Under appropriate consistency and 
stability assumptions it is possible to derive formula ([4]) from this result (as done in 
[31)] ) instead of proving it by the application of the classical Chernoff's Theorem to 
evolution semigroups. It is however amongst our aims to emphasize that semigroup 
techniques may be used to prove approximation results also for non-autonomous 
problems. 

Remark 2.6. In case B(t) = 0, we recover the well-known representation formula 

n-l 

U(t,s)x= lim Ue^^+^^x, 

p=0 
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see Nickel [301 Proposition 3.2] and Schnaubelt [33J Theorem 2.1]. Again, the 
stability condition reduces essentially to the classical stability condition of Kato 

Remark 2.7. It is straightforward to check that if one of the equations is au- 
tonomous, e.g., A(t) = A, then we arrive at the same product formula but we 
can split the original operator C into two (and not three) operators, namely into 
+ A and B(-). 



3. Operator splitting 
In this section we assume that we can solve the non- autonomous equations 
(Eq. A) ±u(t) = A(t)u(t), 

(Eq. B) £„(t) = B(t)v{t) 

and want to construct the solution of (|NCP[) applying an operator splitting pro- 
cedure. For the sake of simplicity we only present the case of sequential splitting: 
We start with the initial value u sq {s) = x, then solve the first equation on the time 
interval [s, s + r]. Then we take this (s + t) as the initial value for the second 
equation which we solve on [s, s + t]. With this result u sq (s + r) := 1*2 (5 + r) as 
initial value for ( |Eq. A[ | we restart the procedure and iterate it n times. Formally: 

J f t u[ k) (t) = A{t)u[ k) (t), te (s+(fc-l)r,s + fcr], 

\ u^(s + (k - l)r) = u s %s + (k- l)r), 

/ = B it)u£\t), te (s+(k-l)T,s + kr], 

\u { 2 k \ s + (k-l)T) = u (k) ( s + kT), 

u sq (s + kT) :=74 fc) (s + fcr), 

for k = 1, 2, . . . , n. If U and V denote the evolution families solving the above 
equations (Eq. A)-(Eq. B), then we have 

uf\r) = U(r, s + (k- l)r) M sq (s + (k - l)r), 

and 

u [k) (r) = V(r, s + (k- 1)t)m^ } (s + fcr) 

= V(r, s + (k- 1)t)U(s + fcr, s + (fc - 1)t)w sc1 (.s + (k— l)r). 

By this the splitting solution it sq can be written as 
fc-i 

u sq (s + fcr) = Y[ (^( s + (P + l ) T > s + pr)U(s + (p+ l)r, s + pr) 

In the following we analyze the convergence of this procedure. 
Assumption 3.1. Suppose that 

a) the non-autonomous Cauchy problems corresponding to the operators A(- 
£>(•), A(-), and B{-) are well-posed, and that 

b) (Stability) there exist M > 1 and w£l such that 



.-eJ 



sup|| f[ Via-frs-te&pia-*^ 

'p — TL—1 



< Me"'. 
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Here, again, the evolution family solving the Cauchy problem corresponding to 
A(-) and B(-), will be denoted by U and V, respectively. Further, we denote the 
evolution family solving (|NCP[) by W and the corresponding evolution semigroup, 
generated by the closure of C = + A(-) + B(-), by W. 



Theorem 3.2. Under Assumvtions \371\ one has the convergence 

ra-l 

W(t, s)x= lim Y[V{s + (p+1 ^- s) , s + &=£p(s + (p+i y~ s) , s + s ^ 1 )x 



n— ¥ oo 

p=0 



for all x £ X . 

Proof. In the space X, we define 

Ht) -=V{;--t) 

and G(t) :=U(-,- -t)K(t). 

Inductively, one can see that 

{Hmi)) n f={n-r-i)u{-,--i)n{i)) n f 

p— n— 1 

By our assumptions, the closure C of the operator C = — -^+A(-)+B(-) generates 
a strongly continuous semigroup on X, and hence the set (A — C)D(C) is dense. 
Straightforward calculation analogous to the one in the proof of Theorem l2.2l vields 
that (T(-)Q(-)) (0)/ = Cf for / G D(C). Hence, by the stability assumption, we can 
apply Chcrnoff 's Theorem to this function and obtain that the evolution semigroup 
generated by C is given by 

p—n—l 

From this, by picking some / £ X with /(s) = x, we obtain for the evolution family 
W(t, s)x = 



lim TT v(t-&=^,t- »+ i y- s ) )c/(t-e^£l,f 

,„vrsn A A \ 71 71 / V 71 



(p+l)(t-s) \ 



n—too 

p—n—l 



= l^]Jv(s + ^^,s + ^)U(s + ^^^ 

p=0 

which was to be proved. □ 

Remark 3.3. Note that the stability condition is trivially satisfied if the evolution 
families U and V are quasi-contractive, i.e., if AI < 1 can be taken in Definition 
11.51 (iii). In general, as usual with stability assumptions, it is rather hard to verify. 

Using similar arguments but a different decomposition, we arrive at a different 
splitting formula using evolution families corresponding to different (timc-rescaled) 
evolution equations. 

Proposition 3.4. Suppose that the operator families A(-/2), B(-/2) and A(-)+B(-) 
generate the evolution families U, V andW , respectively. Assume furthermore that 
there is M > 1 and wel such that 

supll ff V{2s-^,2s-Qs+^)U(2s- ^S±Dl,2s- ^±^) < Me" 4 
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Then we have 
W[t,s)x 



n-l 



= lim TT V(2s + WX'-') , 2s + ^ +1 ^ )U(2s + , 2s + ™^)x. 

Proof. In the space A", we write formally 

+ ^ + 5(0 = ( ~ 2^ + A ^ + (-2^ + = * + Bl - 

Since the division by 2 in the formula means a rescaling of the corresponding evo- 
lution semigroups S and 7~, we obtain the representation formulas 

S{t) = V(2-,2- -t)K{t/2) 

T{t) = U(2-,2- -t)K(t/2). 

By induction one can see that 

{S(i)m)Yf = (V{2; 2 • -i) n{t/2n)U{2; 2 • -±)K(t/2n)) n f 

= Yl V{2 ■ 2 • - (2p + 1)t )tT(2 ■ - (2p + 1}t , 2 • - (2p + 2)t )/(- - t). 

p— n — 1 

Again, the closure C of the operator C = — 4- + A(-) + B(-) generates a strongly 
continuous semigroup on X, hence (A — C)D(C) is dense. By this and by the sta- 
bility assumption Chernoff 's Theorem is applicable. We obtain that the evolution 
semigroup generated by C is given by 

W(f)/= Hm f[ V(2--^,2--£^)U(2-~£^,2--&±?>l)f(.-t). 

p—n—l 

By passing to the evolution family we get the assertion: 
W(t,s)x 

= lim T\y(2t- 2 S^,2t- ^ +1 ^ )U(2t- (2p+D(*-*) 2 f - Vp+W-) \ x 



n— >oo 

p—n — l 

n-l 



n— >oo 

p=0 



lim TT y (2a + , 2 S + ^ +1 ^ )u(2s + WW-') , 2s + ^^)x. 



□ 



Remark 3.5. Note that, in contrast to the autonomous case, there is no general 
connection between the evolution families U and U, see Nickel [29] . 

4. Generalizations and Remarks 

Higher order splitting methods. We now show how the previous results gen- 
eralize to higher order splitting methods. The results are, using the stage set up 
previously, direct applications of the corresponding autonomous results applied to 
the evolution semigroups. We restrict ourselves to the Strang and symmetrically 
weighted splitting, but other splitting methods can be handled analogously. In 
any case only the stability condition has to be adapted. This, however, is always 
satisfied (and typically verifiable) if the operators involved are contractions. 

Theorem 4.1. Suppose that Assumvtions [K7\ a) and c) are satisfied, and that the 
stability condition holds in the following form: 
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< Me" 



b>) sup ff e^-^e^-^e^"^ 

seR " P i„-i 

in the case of the Strang splitting, or: 

b ») SU pi|| A UMs-*-k) e iB { s-?4) +e iB(s-EL) ei A(s-^ 

s6 r2™II 11 A 
p— n— l 

in the case of the symmetrically weighted splitting. Then we have 

n~l 



< Me"* 



W(t,s)x = lim TTe^^+^^e 



'e 

n— Kx} - 1 A 

for all x E X in case of the Strang splitting; and we have 
W{t,s)x 

1 n — 1 

n^oo 2 71 11 V 
p=0 

for all x £ X in case of the symmetrically weighted splitting. 

Proof. The statements follow immediately by the same reasonings as in the proof 
of Theorem 12.21 but now considering the expressions 

for the Strang-splitting, and 

I( K (i)(e^Oe^(-) + e^(-)e^(-)))", 

for the weighted splitting, respectively. □ 

Remark 4.2. It can be shown by exactly the same arguments as in Csomos and 
Nickel [7J Lemma 2.3] that the stability condition (b') is equivalent to the stability 
condition in Assumption 12 . 1 1 b) for the sequential splitting. 

Spatial approximations. Continuing earlier investigations started in Batkai, 
Csomos and Nickel [2], we show that operator splitting combined with spatial ap- 
proximations is also convergent. We only concentrate on the formula ([!]) for the 
sequential splitting. Other methods can be considered analogously. 

Assumption 4.3. Let X m , m £ N be Banach spaces and take operators 

P m : X — > X m and J m : X rn — > X 

fulfilling the following properties: 

(i) P m J m = I m for all m £ N, where I m is the identity operator in X m , 

(ii) lim J m P m x = x for all x G X, 

m— >oo 

(iii) || J m 1 1 < K and ||P m || < K for all m € N and a suitable absolute constant 
K > 1. 

The operators P m together with the spaces X m usually refer to a kind of spatial 
discretization method (triangulation, Galerkin approximation, Fourier coefficients, 
etc.), the spaces X m are in most applications finite dimensional spaces, and the 
operators J m refer to the interpolation method describing how we associate spe- 
cific elements of the function space to the elements of the approximating spaces 
(linear/polynomial/spline interpolation, etc.). 
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Assumption 4.4. For each m G N and r G K let the operators j4 TO (r) and B m (r) 
be generators of strongly continuous semigroups e tAm ^ and e tBm ^ r \ respectively. 
Assume furthermore that 

a) (Stability) there exist constants M > 1 and wet such that 

\\ e hA(r)^ \\ e hA m (r)^ \\ c hB(r)^ || e fcBm(r)|| < M{J «k for ^ ]% > g ^ f £ R) ^ 

supll fr (e^-^e*^-^) < Me"', and that 

b) (Consistency) the identities lim J m A m {-)P m f = A(-)/ for all / G _D(A(-)), 
and lim J m B m {-)P m f = B(-)f for all / G £>(£(•)) hold. 

m— >oc 

As in Batkai, Csomos and Nickel [2], stability and consistency implies conver- 
gence. 

Theorem 4.5. Suppose that Assumption ^. J\ is satisfied. Then one has the con- 
vergence 

n-l 

W(t,s)x= lim limJ^fc^^^'e^^^*)^ 



m— >oo n— >oo 

p=0 



for all x G X . 



Proof. We will apply Batkai, Csomos and Nickel Theorem 3.6], the modified 
Chernoff 's Theorem directly. To this end, define the spaces 

X m = Co(R; X m ), X :=Co(R;A) 

and the projection operators 

V m = I®P m :X^ X m , (V m f)(t) := P m f(t), 

and interpolation operators 

iJ in I ® Jrn • X m Y X , \t!Tmfm)\^) ■ — Jmfm(t*). 

We have to check that these operators satisfy the conditions in Assumption 14.31 
Conditions (i) and (iii) are immediate from the definitions. The (JrriP m f)(s) — > 
f(s) is true pointwise. We have to show that the convergence holds in fact uniformly 
in s G R. Take e > 0. Let / G X and [a, b] C K such that ||/(s)|| < ^2 for all 
s GM \ [a,b]. Then 

\\J m P m f(s)-f(s)\\<e 

for s G K \ [a, b]. Since / is uniformly continuous, there is 5 > such that for 
all s, t G [a, 6], |s — i| < <5, we have \\f(s) — f(t)\\ < K i +2 ■ Take a partition 
a = So < s± < . . . < s n = b such that |sj+i — Si\ < S. Then by definition, there is 
M > such that for all m> M 

\\J m P m f(si) - f(8i)\\ < jfa. 

Since for s G [a, 6] there is j such that s G [sj, Sj+i], we get for m > M, 

||J ro P ro /(s)-/(s)|| 

< \\J m P m (f(s) - /(* 3 -))|| + \\J m P m f( Sj ) - /(* 3 -)|| + \\f(sj) - f(s)\\ < e. 

Hence, \\J m Pmf — f\\oo < £ holds for all m > M . 

The validity of Assumption 14.41 implies that the corresponding multiplication 
semigroups satisfy the necessary stability and consistency conditions. □ 
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Positivity preservation. As it was pointed out by W. Arendt (Ulm), the prod- 
uct and splitting formulas can be used to show positivity properties of evolution 
families. On the terminology and properties of positive operator semigroups see 
Arendt et al. Q] or Engel and Nagel [8j Section VI. 1]. 

Theorem 4.6. Assume that X is a Banach lattice. 

(1) Let the conditions of Assumvtions \2. 1\ are satisfied and that all the operators 
A{r) and B{r) generate positive semigroups. Then the evolution family W 
given by ^ in Theorem \2.2\ is vositive. 

(2) Let the conditions of Assumvtions \S.l\ are satisfied and that all the evolu- 
tion families U and V are positive. Then the evolution family W given by 
Theorem \3.2\ is positive. 



The proof is an immediate consequence of the fact that the corresponding mul- 
tiplication, shift, and evolution semigroups are positive. It would be an important 
and interesting question whether similar results hold for shape preserving semi- 
groups in the sense of Kovacs [HJ Definition 20] . 

5. A NON-AUTONOMOUS PARABOLIC EQUATION 

In order to demonstrate the range of our results, we will consider an important 
and much studied parabolic equation 

(5) d t u{x, t) = Au(x, t) + V(x, t)u(x, t) 

in R d with appropriate initial conditions, where V is a smooth and bounded func- 
tion. Rewritten abstractly this takes the form 

(6) &u(t) = Au(t) + V(t)u(t) 

with u : R+ —> L 2 (R d ) =: X a vector valued function. Hence a straightforward 
choice for the splitting for the evolution semigroups is 

A := — 4- + A, B := the pointwise multiplication by V(t). 

These operators (with appropriate domain) generate the following semigroups on 
the Banach space X := BUC (R; L 2 (R d )) 

\T{t)f]{s) := c tA f(s - t) and [S(t)f](s) := e tv ^f(s). 

We shall assume that V € BUC (R; L°°(R d )), so B is bounded. The domain of 
the generator of S can be given explicitly, see Nagel, Nickel and Romanelli [2"51 
Prop. 4.3]): 

D(A) = {/ G BUC (R; X) n BUC J (R; X_i) : -/' + A_i/ G BUC (R; A)}, 

here A_i with domain L 2 (R d ) is the generator of the extrapolated semigroup, see 
Engel and Nagel [5J Section II. 5. a] for the corresponding definitions. 

As a corollary of Theorem 12.21 we obtain the convergence of the sequential (and 
also the Strang) splitting procedures. 

Proposition 5.1. Suppose that the potential V E BUC (R; L°°(R d )). Let W denote 
the semigroup generated by A + B on BUC (R; L 2 (R d )). For every function f G 
BUC (R; L 2 (R d )) we have the product formula 

hm (5(1)T(£)) B / = W(t)f, 

where the convergence is uniform on compact time-intervals. Let (W(t, s)) t > s de- 
note the evolution system solving ([6]) on L 2 (R d ). Then for every uq G L 2 (R d ) we 
have 



lim 



W(t, s)u - Y[ e i ^ £y(s+ ^ ) e^ £A M 

p=0 



0, 
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locally uniformly for s <t. 

Proof. For the first assertion we only have to verify the stability Assumption 12.1 
b), and then the assertion follows directly from Chernoff's Theorem ll.il Stability 
follows, because the semigroup e tA is contractive and V(s) is uniformly bounded. 
The second assertion is a direct consequence of Theorem 12.21 □ 



Next we study convergence rates for the sequential splitting procedure applied 
to the above equation ([5]). To this end we apply Theorem II .21 to the corresponding 
evolution semigroups. 

Theorem 5.2. Suppose that V G BUC (R; W 2 '°°(R d )) n BUC J (R; L°°(R d )). If 
f e BUC^R;!! 2 ^)), we obtain 

Ct 2 

ll(S(£)T(£)) n -W(t)/|| < _||/|| BUC1(R;H2(Rd)) . 

Before we prove the theorem, let us first reformulate this product formula for 
the solutions of the non-autonomous problem. 

Corollary 5.3. Consider the non- autonomous parabolic equation 

jd t u(x, t) = Au(x, t) + V(x, t)u(x, t), t>s, x e R d , 
}u(x,s)=u (x), x € R d . 

Suppose that V <E BUC (R; W 2 >°°(R d )) H BUC L°°(R d )). If u G H 2 (R d ) then 
for the evolution family (W(t, s))t>« solving the above problem we have 

n-l 

UQ 

p=0 



IU(t, s ) uo -II e ^ y(s+ " )e ^ 



< N WW- 



Proof. The assertion follows from Theorem 15.21 from the calculations in the proof 
of Theorem 15.11 and from the fact that the constant function f(s) := uq £ H 2 (R d ) 
is in the domain of A. □ 



In order to prove Theorem 15.21 we have to verify the commutator condition in 
Theorem 11.21 for the generators of the evolution semigroups. To do this, we need 
the following abstract identification of the domains of fractional powers of evolution 
semigroup generators. 

In what follows, let B (R; Y), BUC "(R; Y) etc. denote the space of bounded Y- 
valued functions, the space of a-H61der continuous Y- valued functions etc., where Y 
is some Banach space. Let X be a fixed Banach space, and let e be a (contractive) 
analytic semigroup with generator (A, D(A)) thereon. The fractional powers of —A 
are denoted by ((— A) a , D((— A) Q )). Denote by F a the abstract Favard spaces for 
X and (c tA )i>o, i.e., 

F a := \ x G X : ||x|| Q := ||a;|| + sup II ° tA l~ x II < +oo|, 

which becomes a Banach space if endowed with the norm || ■ || a . For every a, /3 G 
(0, 1) with a > (3 we have continuous embeddings (see Engcl and Nagel [3 Sec. II. 5.]): 

F a ^ D((-Af) ^ Fp. 

Consider now the Banach space X := BUC (R; X) and the semigroup 

(T(t)f)(s) :=e tA f(s-t) 
thereon. We are interested in the Favard spaces X a of this semigroup. 
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Proposition 5.4. In the above setting we have the following continuous inclusions: 
BUC (R; D((~A) a )) nX a ^ BUC (R; D((~A) f3 )) n BUC^R; X), 

/or all < /3 < a < 1, and 

BUC Q (M; X) n BUC (K; L>((-A) Q )) c — > BUC (R; L>((-A) /3 )) n Xp, 

for all < j8 < a < 1. 

Proof. We show the statement for j3 = a, the rest then immediately follows. We 
start with the second inclusion. For / G BUC (R;X) we can write 



sup 

t>o 



?W - / l|c tA /(« - *) - fjs) 
= sup sup 

t>0 s6K <Q 



fc° 



_ mm mm l|e M /(^)-/( S ) + c^(/( g -t)~ /( 5 ))|| 

— S LI JJ o LI J_) 

t>0 sGK 

< SUp ||/(s)||F a + II/HbUC"- 

This shows that if / G B (R; F a ) n BUC Q (R; X), then / G and the inclusion 
is continuous, i.e. 



\\f\\ Xa <c(\\fh 



l;F a ) ~r \\J \\BVC a (R;X] 



To see the first inclusion we use now that A generates an analytic semigroup. If 
/ G BUC(M;L>((-.4) Q )), then 

|jc^/(^t)-/(.-t)|| _ \\(c tA -I)(-A)-«(-Arf( s -t)\\ 

□ UU — - o LL L) 

t>o t a t>o t a 

<Csup||(-A) Q /( s -<)|| <C\\f\\ BUC(R;D((-A)»))- 

This implies then 

sup imi^kjm < supsup iM-i +cn/ii BUC(R;B(( _ A)(l)) . 

t>o t seK t >o t 

The proof is complete. □ 

Now we are in the position to check the required commutator condition and thus 
to prove Theorem[ 



Proof of Theorem \5.2i Consider now the evolution semigroup corresponding to the 
non-autonomous equation (J5j> . The corresponding generator is given formally as 

-& + A + V(t). 

Take now / G BUC 1 (R;H 2 (R d )), and notice that then / belongs to the domain 
D(A). We calculate the commutator of A and B. We have 

[A,B]f = -V'{t)f{t) + (AV(t))f + 2W(t) • V/(t). 

Now, if we assume that V G BUC ^R; L°°(R d )) and V G BUC (R; W 2 >°°(R d )), then 
the first two terms can be estimated by c||/||, so we have only to deal with the term 
2VU • V/, for which it suffices to estimate dif(t) for i = 1, . . . , d. We have 

\\dif(t)\\ 2 < c\\A l ' 2 f(t)\\ 3 (di is A^-bounded on L 2 ). 

By Proposition 15.41 this completes the proof of the commutator condition (TTJ in the 
form 

< ||(-^) Q /|| for all / G D(A) with some given a > 1/2. 
Hence Theorem 11.21 yields the assertion. □ 
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6. Numerical examples illustrating the convergence 

In Section [5] we already introduced the non- autonomous parabolic equation 
(sometimes also called imaginary time Schrodingcr equation) 

d t u(x, t) = Au(x, t) + V(x, t)u(x, t) 

in M. d with appropriate initial conditions with V being a smooth and bounded func- 
tion. In the following we will apply the sequential splitting introduced in Section [3] 
to the sub-operators 

A(t) := A and B(t) := multiplication by V(x,t). 

In Theorem 12.21 we showed that the product formula describing the sequential 
splitting is convergent also in the case if we are able to solve the corresponding 
autonomous Cauchy problems (Eq. l)-(Eq. 2) with operators A(r) and B(r) for 
every time level r £ R. We will use this result when constructing our numerical 
scheme. 

In order to illustrate numerically the convergence of the sequential splitting and 
give an estimate on its order, let us consider the following non-autonomous equation 
with boundary and initial conditions: 

t) = d 2 x u{x, t) + V{x, t)u(x, t), t > 0, x e [0, 1], 
(7) iu{0,t)=u(l,t) = 0, t>0, 

u(x,0) = u (x), x e [0, 1] 

with functions V(x,t) and uq(x) given later on in the example. 

6.1. Error analysis. Let (u sp i)™ denote the approximation of the exact solution 
u(iS, nr) of problem (JT]) at time nr and at the grid point i5 (with n = 0, N — 1 
and i = 0,...,I — 1) using sequential splitting. At this point the time-step r = 
j^—f and the grid size S = have certain given values. We call (ix sp i)" = 
((Wspi)oi ( u s P i)i, (w sp i)™_i), n = 0,1, ..,N — 1. the split solution of problem 0. 
As already seen, the order of the splitting procedure can be estimated with the help 
of the splitting error defined by 

:= K-< pl || 

where u n = (uq, u", u™^) with uf = u(i6,nr), i = 0,1,...,/— 1. With this 
notation the splitting procedure (or an arbitrary finite difference method) is of 
order p > if for sufficiently smooth initial values there is a constant C > such 
that for all t G [0, to] we have 

(■n < c 

c spl ^ n P I 

or, if the method is stable, equivalently, 

In general, the exact solution of problem ([7]) is unknown, therefore, the local split- 
ting error j is to be estimated as well. To this end we compute a so-called 
reference solution w™ ef on a finer space grid using no splitting procedure. Then 
the order p of the splitting procedure can be determined as follows. From the 
definition of p we have £^ pl < Ct p+1 . Approximating u n with u" f , we obtain 

#pi « £s P i : = Kef - <ill < Cr P +\ Thus, 

log§ p i < (p+l)logr + logC7. 

Then we can estimate p by computing the approximate local splitting error £ spl for 
many different values of the time-step r, plotting the logarithm of the results, and 
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fitting a line of form y(w) = aw + b to them. Hence, a ss p + 1 and b « logC. 
Note, however, that the split solution contains not only the splitting error but also 
a certain amount of error originating from the spatial and temporal discretization. 
In what follows we show how to determine the numerical solutions wj ef and u^pi- 

We also note that it is reasonable to compute a relative local error defined as 

f 1 

F - SP 1 

<-"loC ,i 1 |i 

Kef II 

because this yields the ratio how the split solution differs from the reference solution. 

6.2. Numerical scheme. In order to solve numerically the problem Q we should 
discretize it in both space and time. For the temporal discretization we used the 
Crank-Nicholson method, and we chose the finite difference method for the spatial 
discretization. 



6.2.1. Reference solution. As mentioned above, we need a reference solution u" ei 
computed without using splitting procedures. After discretizing the equation, we 
obtain the following numerical scheme for determining (u re f)™ +1 : 

(8) Wef), = (1 - (i/ r ef)r +1 ) _1 (l + (# r ef)r)«ef)i 

with 

t ( v n+1 - 2v n+1 + 7/" +1 \ 

l-ttrefjj - g I p + V i U 

where V™ := V(iS, nr). 



6.2.2. Split solution. Application of sequential splitting means that instead of the 
whole problem (J7]) two sub-problems are solved. In our examples the first sub- 
problem corresponds to the diffusion equation d t UA{x,t) = d%UA(x,t). Its numer- 
ical solution u n A can also be computed using Crank-Nicholson temporal and finite 
difference spatial discretization methods. Then we obtain the following numerical 
scheme similar to flSJ): 

(9) = (1 - (H A )- +1 y\l + (H A )?)(u n A )i 



with 

* - 2 P ' 



The second sub-problem has the multiplication operator by V(x, t) on its right-hand 
side, i.e. dtus(x,t) = V(x,t)uB(x,t). We refer again to Theorem 12.21 and take the 
function V only at time levels t = nr, n = 0, 1, N — 1. In this (autonomous) 
case the exact solution us(x,t) = e tv ^ x ' nT ^u^{x) is known. At the n th time level 
and on the space grid it has the form 



(10) 



(u n B )i = u B (i6,nT) = e rV ^ n ^u (id). 
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Due to the product formula ©, the split solution u™ pl is given by the following 
algorithm: 

for i = Q,. ..,1-1 

initial function: (ua); := uo(iS) 

end 

for n = 0,1,..., N- 1 
for i = 0,1,..., I- 1 

solve the first sub-problem using © ==> (it a)™ 

end 

for i = 0,1,..., I- 1 

solve the second sub-problem using (fT0|) => (ub)? 

end 

end 

split solution: u^f 1 := Mg -1 

6.3. Numerical results. Now we present some numerical results on the following 
example. 
Choose 

V(x,t) = t- 500a; 2 and u (x) = e" 50 ^" - 4 ) 2 . 

Since the exact solution is unknown in this case, we should estimate the local 
splitting error using the reference solution instead of the exact one. Then the 
relative local splitting error £\ oc and its order p can be measured. 




FIGURE 1. Numerical solution of equation 10 at time levels t = 0, 
t = 10 -3 , t = 5 ■ 10~ 3 , and t = 10 -2 , respectively. 



On Figure [T] the time-behavior of the reference solution can be seen at the four 
time levels t = 0, t = 10~ 3 , t = 5 • 10~ 3 , and t — 10~ 2 , respectively. The effect 
of the diffusion can be clearly observed. Figure [5] shows the result of the fitting. 
The dots correspond to log(£i oc ) for the various step sizes. The line fitted to these 
points has the form ?/(log(r)) = alog(r) + b with a = 1.9470 and b = 3.25925. 
As mentioned above, the order of the splitting procedure p can be estimated by 
a — 1 « 1, that is, the sequential splitting is of first order. 
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-4 -3.8 -3.6 -3.4 -3.2 -3 -2.8 -2.6 -2.4 -2.2 -2 
log(T) 



FIGURE 2. Results obtained by applying the sequential splitting with 
various time steps (dots), and the line y(w) = aw + b fitted to them with 
parameters a = 1.9470 « p + 1 and b = 3.25925. 
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